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It is often stated that gamma-ray bursts (GRBs) have typical en- 
ergies of several hundred keV, where the typical energy may be char- 
acterized by the hardness H, the photon energy corresponding to the 
. . . peak of vFv. Among the 54 BATSE bursts analyzed by Band et. al. 

more than half have lOOkeV < H < 400keV. Is the narrow range of 
H a real feature of GRBs or is it due to an observational bias? We 
consider the possibility that bursts of a given bolometric luminosity 
occur with a distribution: p{H)dlog H oc H~'d\ogH. We model the 
O , detection efficiency of BATSE as a function of H and calculate the ex- 

pected distribution of H in the observed sample for various values of 7. 
' The Band sample shows a paucity of soft (X-ray) bursts, which may be 

real. However, because the detection efficiency of BATSE falls steeply 
with increasing H, the paucity of hard bursts need not be real. We 
find that the observed sample is consistent with a distribution above 
H — 100 keV with 7 ~ (constant numbers of GRBs per decade of 
^ ' hardness) or even 7 = 0.5 (increasing numbers with increasing hard- 

ness). Thus, we suggest that a large population of unobserved hard 
(N . gamma-ray bursts may exist. It is important to extend the present 

analysis to a larger sample of BATSE bursts and to include the OSSE 
and COMPTEL limits. If the full sample is consistent with 7 > 0, then 
it would be interesting to look for MeV bursts in the future. 
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One striking feature that is common to all gamma-ray bursts (GRBs) is the 
fact that most of the observed photons correspond to low energy gamma-rays, 
with energies of a few tens to few hundreds of keV. While other features of the 
^ bursts, in particular the temporal structure, vary significantly from one burst 
, to another, this feature seems to be quite invariant. One wonders, therefore, 
^ • whether this is a clue to the nature of GRBs — a phenomenon that theorists 
• i-H , should strive to explain — or if it is just the consequence of an observational 
' bias against detection of harder or softer bursts. 

Band et al. (|l|) have analyzed 54 strong GRBs (hereafter the Band sample), 
fitting the spectra using a four parameter function: 
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k L lookcv J ^MP-a)iimcv) ^ ior E > {a - f3)Eo. 
This function, which provides a good fit to most of the observed spectra. 
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is characterized by two power laws joined smoothly at a break energy H = 
(a — (3)Eq. For most observed values of a and (3, vF^ cx E'^N{E) rises below 
the break and decreases above it. The break energy H is thus the "typical" 
energy of the observed radiation, in the sense that this is where the source 
emits the bulk of its luminosity. H is correlated with, but not equal to, the 
hardness ratio which is commonly used in analyzing BATSE GRBs, namely 
the ratio of photons observed in channel 3 to those observed in channel 2. 

Figure 1 shows the distribution of observed values of H in the Band sample. 
There is a clear and marked maximum in the distribution for H ^ 200keV, 
and most of the bursts lie over the range 100 keV < H < 400 keV. Should we, 
therefore, conclude that most GRBs have hardnesses around 100-400 keV? 

To answer this question we have used a simple model of the BATSE detector 
to calculate the expected hardness distribution of GRBs detected by BATSE 
for various assumed intrinsic hardness distributions. We have then compared 
the expected and observed distributions to see which intrinsic distributions 
are consistent with the data and which are not. 

We assume that the intrinsic hardness distribution is as follows: 

• All GRBs are adequately described by a Band spectrum, characterized 
by the parameters, a, (3, and H. 

• The number of GRBs in a logarithmic interval, dlogH, varies as: 
p{H)d log H = poH^d log H. 

• We also allow for the possibility that the intrinsic luminosity, L, is cor- 
related with H: LoiH/Ho)^. 

We now calculate, for given 7 and ^, the distribution of H values of bursts 
seen by a detector like BATSE. We make use of the fact that BATSE triggers 
on counts in the second and third energy channels, with photon energies 
between 50keV and 300keV. We calculate, therefore, how many bursts with a 
given H yield a count rate in the 50keV to BOOkeV band that is larger than a 
threshold, Cmm- (We ignore the fact that BATSE's counts in these channels 
may correspond to higher energy photons and that the intrinsic spectrum is 
related to the observed counts via the DRM, (||)). The count rate depends 
on a and /3, in addition to the above mentioned dependence on ^ and 7. We 
use the average values given by Band et al. (0): a = —0.73 and /3 — —2.22. 
For reasons that will become clear shortly we also use a — —0.41, which is 
the average value of a for GRBs with 50keV < H < 300keV. 

The count rate in the 50-300keV energy band from a burst at a distance D 

is: 
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where A is the area of the detector. The quantity C represents the number of 
counts for every erg of energy incident on the detector (counts/ergs), and is 
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easily calculated once the shape of the spectrum {a, (3, H) and the detector 
limits (50, 300 keV) are given. The second equality follows from the assumed 
correlation between L and H. For simplicity we use here and in the following 
a Newtonian geometry and ignore cosmological effects. This is justified since 
the Band sample is composed mostly of strong bursts for which cosmological 
redshift and spatial curvature effects are small. A detailed discussion that 
includes these effects will be published elsewhere (Cohen, Narayan & Piran, 



The number of bursts with hardness H detected by a detector with a lim- 
iting sensitivity Cmin is: 



Figure 1 depicts N(H). the expected munbcr of detected bursts in a loga- 
rithmic interval of H, as a function of \og{H). We have shown the calculated 
curves for a = a and a = a and for 7 + 3^/2 = and 7 + 3^/2 = 0.5; in all 
cases, we set (3 = 0. The observed distribution corresponding to the Band 
sample, with statistical error bars, is also shown on the same plot. 

A comparison of the four calculated (or expected) curves with the observed 
distribution reveals immediately a paucity of soft bursts. In all four cases, the 
number of soft bursts we would have expected to see is significantly larger than 
the number actually observed. Therefore, unless BATSE has an unexpectedly 
large selection bias against detecting soft photons (i.e. significantly poorer 
sensitivity than we have assumed in our model), we conclude that the lower 
cut-off in the observed distribution of hardnesses is a real phenomenon. That 
is, there really are very few soft GRBs, and N{H) does have a lower cutoff. 

The story is, however, very different for larger values of H. The data show 
very small numbers of hard bursts; for instance, only two bursts out of 54 have 
H > IMeV. Nevertheless, this docs not mean that GRBs intrinsically cut-off' 
for hardnesses above IMeV. For instance, our theoretical model with ^ = 
and 7 = 0, with equal number of bursts per logarithmic hardness interval, 
actually predicts that even fewer bursts should have been observed. The best 
fit to the data is obtained with 7 + 3^/2 = 0.5 which corresponds to a burst 
population with an increasing number of bursts with hardness (for ^ = 0) or 
to an equal number of bursts per logarithmic hardness interval (7 = 0) but 
with harder bursts being more luminous (^ = 1/3). 

The interpretation of the result is quite simple. There is an observational 
bias against detecting bursts with H > 500keV by current detectors. Two 
factors operate. For bursts with a fixed luminosity, harder bursts have fewer 
photons. This makes the detection of harder bursts more difficult in any 
detector that is triggered by photon counts. The decrease in sensitivity in 
BATSE is even more severe since BATSE triggers on photons in the 50keV to 
SOOkeV range and as the bursts becomes harder most of the emitted photons 
are further and further away from this energy range. 
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FIG. 1. Observed and expected numbers of detected bursts per interval of \og{H). 
The the error bars represent statistical errors. The best fit (solid line) corresponds 
to a = a = -0.41, 7 + 3^/2 = 0.5. Other curves are: a = a = -0.7, 7 + 3^/2 = 0.5 
(dotted line), a = a = —0.41, 7 + 3{/2 = (dashed line), and a = d = —0.7, 
7 + 3^/2 = (dashed-dotted line). 



The last point depends, of course, on the power law index a in the low 
energy range. It suggests that hard bursts that are detected by BATSE will 
have values of a lower than average. Fig. 2 depicts the distribution of a 
values in different hardness regimes, and shows the effect clearly. Among the 
54 bursts in the Band sample, hard bursts do have significantly more negative 
values of a. It is for this reason that we have introduced a, the average a value 
over the intermediate hardness range, 50 keV < H < 300 keV. If the intrinsic 
a distribution is independent of the hardness than a is a better estimate of 
the average a than a. 

Figure 2, by itself, without any of the theoretical arguments presented ear- 
her, suggests the existence of the selection effect we have discussed. There 
are two ways of interpreting the evidence in this plot. One could say that the 
plot reflects the true distribution of burst properties and that for some reason 
a happens to be correlated with H in the particular manner seen in the data. 
This is very ad hoc. The alternative, which we find much more attractive, is 
to say that BATSE has difficulty detecting hard bursts, and finds it particu- 
larly difficult to detect hard bursts with less negative, or positive, values of a; 
that is, such bursts do exist, but BATSE misses them. If we accept the latter 
explanation, then it means that BATSE is certainly missing at least some part 
of the hard population, namely those bursts with less negative (or positive) 
a. It is then but a small step to accept the entire argument presented earlier. 

Our main result then is the following. When observational selection effects 
are taken into account, the two hard bursts with H > IMeV in the Band 
sample are but "the tip of the iceberg," and represent a large number of 
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FIG. 2. The low energy spectral slope a vs. the hardness H in the Band sample. 
Hard bursts tend to have more negative a, which we interpret as evidence for the 
selection effect discussed in the paper. The vertical lines mark 50keV and SOOkeV 

undetected hard GRBs. The observed distribution of hardness is, in fact, 
consistent with the possibility that most GRBs are harder than the bursts 
detected by BATSE, and that there is a large population of mostly undetected 
GRBs whose typical photon energy is in the MeV, or even harder, range. 

So far we have compared the theoretical predictions to the hardness distri- 
bution in the Band sample. It will be interesting to perform the same analysis 
on a larger sub-sample (see (|3|)). It will also be of great interest to perform a 
similar analysis on data from other GRB detectors. One intriguing possibility 
is to search in BATSE's raw channel 4 data for untriggered events, which 
might be some of the missing hard bursts. This would be a complementary 
search to the soft GRB search reported by Kommers et. al. In the mean- 
time, before we know whether hard bursts exist or not, we should be very 
cautious about performing correlations between various characteristics of the 
observed bursts and their hardness parameters. The data are biased, as far as 
H is concerned, and any correlation seen (e.g. the correlation between a and 
H in Fig. 2) might reflect nothing more than selection effects in the data. 
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